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Summary. — We study the formation and detection at the next linear e"'"e~ 
colhder of bound states of level-1 quark Kaluza-Klein excitations Bkk within a 
scenario of universal extra-dimensions (UED). 

PACS 12. 60. -i - Models Beyond the Standard Model. 
PACS 11 . 10. St - Bound states effects. 



1. Introduction 

It is well known that as early as 1921 Theodore Kaluza proposed a theory that was 
intended to unify gravity and electromagnetism by considering a space-time with one 
extra space-like dimension [1]. A few years later Oscar Klein proposed that the extra 
space dimension (the fifth dimension) is in reality compactified around a circle of very 
small radius [2]. These revolutionary ideas have thereafter been ignored for quite some 
time. However recent developments in the field of string theory have suggested again the 
possibility that the number of space time dimensions is actually different from Z? = 4 
(indeed string theory models require D = 11, i.e. seven additional dimensions). In 
1990 it was realized [3] that string theory motivates scenarios in which the size of the 
extra dimensions could be as large as i? 10~^^ cm (corresponding roughly to electro- 
weak energy scale (« TeV) contrary to naive expectations which relate them to a scale 
of the order of the Planck length Lp k, 10"'^^ cm (corresponding to the Planck mass 
Mp = .JtKjG « IQi^ GeV). See also [4]. 

Subsequently two approaches have been developed to discuss the observable effects of 
these, as yet, hypothetical extra dimensions. One possibility is to assume that the extra 
space-like dimensions are flat and compactified to a "small" radius. This is the so called 
ADD model [5] where only the gravitational interaction is assumed to propagate in the 
extra-dimension. A second possibility is contemplated in the Randall- Sundrum type of 
models where the extra dimensions do have curvature and are embedded in a warped 
geometry [6, 7]. 

Universal extra-dimensional models were introduced in ref. [8] and are characterized, 
as opposed to the ADD model, by the fact that all particles of the Standard Model (SM) 
are allowed to propagate in the (flat) extra space dimensions, the so called hulk. Here 
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to each SM particle X^^^ corresponds in this model a tower of Kaluza-Klein states X^") 
(KK-excitations), whose masses are related to the size of the compact extra dimension 
introduced and the mass of the SM particle via the relation m'^,^^ « ^ / ■ 

An important aspect of the UED model is that it provides a viable candidate to the 
Cold Dark Matter. This would be the lightest KK particle (LKP) which typically is 
the level 1 photon. Many aspects of the phenomenology of these KK excitations have 
been discussed in the literature. For reviews see ref. [9, 10, 11, 12]. In particular KK 
production has been considered both at the Cern large hadron collider (LHC) and at the 
next linear collider (ILC). Direct searches of KK level excitations at collider experiments 
give a current bound on the scale of the extra-dimension of the order > 300 GeV. See 
for example ref. [13] . At the Fermilab Tevatron it will be possible to test compactification 
scales up to ^ 500 GeV at least within some particular scenario [14, 15, 16]. 

Lower bounds on the compactification radius arise also from analysis of electro-weak 
precision measurements performed at the Z pole (LEP II). An important feature of these 
type of constraints is their dependence on the Higgs mass. A recent refined analysis [17] 
taking into account sub-leading contributions from the new physics as well as two-loop 
corrections to the standard model p parameter finds that > 600 GeV for a light 
Higgs mass {m,H = 115 GeV) and a top quark mass m/ = 173 GeV at 90% confidence 
level (C.L.). Only assuming a larger value of the Higgs mass the bound is considerably 
weakened down to R^^ > 300 GeV for mn = 600 GeV, thus keeping the model within 
the reach of the Tevatron run II. The finding of this precision analysis are in qualitative 
agreement with previous results [18], but are at variance with the conclusions of a recent 
paper [19] where an analysis of LEP data including data from above the Z pole and two 
loop electro-weak corrections to the Ap parameter pointed to R~^ > 800 (at 95% C.L.). 

It has been shown in ref. [20] that a refined analysis of B ^ Xgj including in addi- 
tion to the leading order contribution from the extra-dimensional KK states, the known 
next-to-next-to-leading order correction in the Standard model (SM) gives a lower bound 
on the compactification radius > 600 GeV at 95% confidence level (CL) and inde- 
pendent of the Higgs mass. 

We discuss here the formation, production and possible detection of bound states 
of Kaluza-Klein n = 1 excitations at e+e" collisions. To estimate the bound state 
contribution to the threshold cross-section, an effect which can be as large as roughly 
a factor of three for strongly interacting particles, we use the method of the Green 
function as opposed to previous works [21] which use a Breit-Wigner approximation. 
The interactions responsible for the formation of level- 1 KK bound states is assumed to 
be an driven Coulomb potential. This allows the use of analytic expressions for the 
Green function of the Coulomb problem. This method has also been recently used by 
the present authors in a study of sleptonium bound states within a slepton co-next to 
lightest supersymmetric particle (slepton co-NLSP) scenario of gauge mediated symmetry 
breaking (GMSB) [22]. 

2. — uiUi Bound State Formation and Production Cross-Section 

In this section we shall review the possible creation of a bound state of the level- 1 KK- 
excitation of the ?i-quark, i.e. a bound state uiUi. The interaction among two Kaluza- 
Klein excitations are driven by the QCD interaction, thus bearing no differences with 
respect to the Standard Model; the strength of the interaction is given by computed 
at a suitable scale [25, 26, 27]. We shall adopt the same formation criterion stated there, 
namely that the formation occurs only if the level splitting depending upon the relevant 
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R-^ (GcV) 


KK mass (GeV) 




Mb (GcV) 


Sis (GeV) 


AE{2P - IS) (GeV) 


400 


478.05 


0.131 


952.57 


3.627 


2.720 


600 


717.08 


0.124 


1429.30 


4.903 


3.677 


800 


956.11 


0.120 


1906.11 


6.089 


4.567 


1000 


1195.14 


0.116 


2382.98 


7.214 


5.411 



Table I. - Results of Coulombic model for the bound state of the level-1 iso-douhlet Ui quark. 
The strong coupling as is computed at the scale Q = r^^ , where rs = 3/{2mas) is the Bohr's ra- 
dius. For each mass value m the scale Q = r~^ depending itself on as must be solved numerically 
from the equation Q = (2/3) mas{Q). 



interaction existing among constituent particles is larger than the natural width of the 
would-be bound state. This translates into the formation requirement 

(1) AE2P-1S > r 

where AE2P-1S = E2P — Eis and F is the width of the would-be bound state. The 
latter is twice the width of the single KK quark, T = 2rKK, as each KK quark could 
decay in a manner independent from the other. 

In our model V{x) is given by a Coulombic potential V{r) = — 4as/(3r) with r = |a;|, 
and where as is the usual QCD coupling constant which has been taken at a suitable 
scale as described in [25, 26]. We are thus able to compute its energy levels given by the 
expression 

(2) £ 

and the separation of the first two energy levels is given by A£'2p_is = ma'^/3. The 

scale at which as is evaluated is given by the inverse of Bohr's radius tb = 3/(2m,as), 
the average distance of the constituents of the bound state. The mass of the nth bound 
state is given by the expression M„ = 2m + £„ where m is the mass of the constituent 
Ml quark and En is given by (2). The wavefunction at the origin, which will be needed 
in order to compute decay widths, for this particular model is given by the expression 
1^(0) |2 = {2mas/3)^/Tv. The results are given in Table I. 

In order to determine whether the bound state will be formed we shall apply the 
criterion given in eq. (1). The if -quark decay widths have been already computed 
in [21], where it has been shown that their values are at most of the order of 100 MeV, 
one order of magnitude less than the energy splittings. In this scenario the eq. (1) 
requirement is always fulfilled, and the bound state is formed for K K-qaaik masses in 
this investigation range. 

In order to describe the cross section of a KK boimd state in the threshold region 
we shall use the method of the Green function. We briefly review here the essential 
features of the mechanism, and refer the reader to the literature for further details [28]. 
Let Qis{x, y, E) be the Green function of the Schrodinger equation which describes the 
bound state by means of a suitable potential V{x). 

The complete expression for the IS Green function of our problem as a function of 
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energy from threshold is given with a shght change of notation by [30] : 



(3) gis{o,o,E + ir) = — 



where k = ^ —m(E + zF). A = 2as?7i/3 and the wave number is = A/A:; E = ^/s — 2m. 
With the position E E + iT we take into account the finite width of the state. The 
ip is the logarithmic derivative of Euler's Gamma function r(a;), 7 ^ 0.57721 is Euler's 
constant and /x is an auxiliary parameter. 

The final expression for the production cross-section of a KK bound state is thus 
given by 

IStt 

(4) a{m, E, T, as) = "^crs Im [Gis] 

where as is the Born expression of the cross-section [31] The process e+e~ U\U\ 
proceeds through the annihilation into the standard model (level-0) gauge bosons 7 and 
Z but in principle one should also consider the contribution of the level- 2 gauge bosons 
7(2) and •^(2)- Especially so in our case of threshold production of the pair u\U\. Indeed 
in this case m « 1/i?, and \/i = 2in + E « 2/R+E and since m^^ ~ 2/i? when producing 
at threshold the UiUi pair we would be close to the 72 and Z2 resonances. However the 
mass spectrum is modified by the radiative corrections. We have verified that over the 
region of parameter space 300 GeV < R'^ < 1000 GeV and 2 < Ai? < 70 the pair 
production threshold 2m„^ is always larger than m^^ 1 ^^22 ^^'^ thus these resonances 
should in principle be included in the calculation. We have also verified, cross checking 
our calculation with the output of a CalcHEP [32, 29] session, that the numerical impact 
of these diagrams is completely negligible. Their contribution turns out to be five orders 
of magnitude; smaller than that of the SM gauge bosons 7, Z. The analytic formula of the 
Born pair production cross section e+e~ — ^ 7*, UiUi can be deduced for example 

from those of heavy quark {tt) [33] taking into account the fact the the level-1 KK quarks 
are vector-like i.e. their coupling to the Z is of the 7^ type and has no axial component. 
Details of the calculation and the explicit expression of the cross section are in [34] . 

In this work we shall concentrate on the continuum region of the cross section, namely 
E > 0. The region below threshold, E < 0, has been already discussed in ref. [21], using 
a Breit-Wigncr description of both the positions and the widths of the peaks. In this 
respect the Green function approach is not expected to point to substantial differences 
relative to the Breit-Wigner one. 

In Fig. 1 we show the cross section for selected values of the scale of the extra di- 
mension, = 400 — 1000 GeV. The results are less sensitive to the other parameter 
(Ai?) which only enters through the logarithmic factors in the radiative correction terms 
in the mass spectrum of the model [34]. In fig. 1 we have fixed Ai? = 20 and varied R~^ 
computing the corresponding values of the level-1 KK quark mass, and assuming the 
energy of the collider being fixed at ^/s = 2mui + E, E being the energy offset from the 
threshold. We have used a value of F = 0.5 GeV for illustrative purpose, compatible with 
the formation of bound state. Different choices of F by even two orders of magnitude 
smaller will not make a visible difference on the figures. 
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Fig. 1. - Production cross-sections of level- 1 KK doublet quark bound states UiUi as a function of 
the energy offset from threshold (y^ — 2mui +E), for values of the scale of the extra-dimension 
i?"^ in the range [400 1000] GeV and a total width of F = 0.5 GeV. The continuous line is 
the Green Function result, the dotted one is the Born approximation given by our analytical 
formula (see [34]). The full circles show the complete agreement with the Born cross section 
from the CalcHEP [29] numerical session including also the annihilation diagrams of 72 and Z2 
whose contribution is however completely negUgible. The cut-off scale A, at which perturbative 
expansions break down, has been fixed so that A.R = 20. 



3. — uiUi Decay Widths 

The KK bound states we discuss here are the pseudoscalar and the vector one 
^Si. For the pseudoscalar state the decay channels are into two photons or two gluons 
for which the following Born level expressions hold (see for instance [26]): 

/c^ T. /lo ^ 4 2487rl^(0)l2 ^ ^ ,1 „ , 2 327rlV'(0)|2 

(5) TbCSo ^ 77) = qta^ and TbCSq ^ gg) = al ^'jj^y' . 

Here qi is the charge of the constituent quark of the bound state, while M and ]^/;(0))^ 
are given by the ones from Coulombic potential previously met. 

The QCD radiative correction results [35], which are the same in the two decays, are 
shown in Fig. 2 (left panel). 

For the vector case the relevant decay channels are the one in charged pairs and 
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Fig. 2. - Left panel Solid line: decay width (kcV) of the pscudoscalar uiui bound state to 
two photons as a function of KK excitation mass. Dashed line: decay width (MeV) of the 
pseudoscalar wiwi bound state to two gluons as a function of the KK mass. Right Panel: 
decay widths of the uiUi vector bound state to two charged particles and into three gluons as 
a function of KK mass. Here we have considered all possible e.m. decay channels. 



the one into three gluons, for which one has 

(6) 

TbCS, ^ q+q^ ) = qUja' and TbC^^i ^ ggg) = ^--^a3__J_^_ 

The charge of the final state charged particle is given by qf. The QCD radiative correc- 
tions [35] depend as usual on , that has to be computed at a scale of the order of 2m. 
The two decays of the vector state are shown together in Fig. 2 (right panel). We observe 
that only the pseudoscalar hadronic decay is in the MeV range and raises approximately 
linearly with KK mass. The ^Sq photonic decay and ^Si decays are smaller by almost 
two orders of magnitude for the considered KK mass range. For the pseudoscalar case 
the hadronic is the dominant decay by far, while in the vector case the decay into charged 
particles, when taking into account all possible processes as seen in Fig. 2 overtakes the 
hadronic decays. 

Other electro- weak decay channels are negligible. Those are proportional to a^, thus 
their ratio to gluonic decays is suppressed by {a/as)"^, at least by two orders of magnitude. 

For most scenarios depending upon the values of A and R [21] single quark decay 
becomes the dominant decay channel for the bound state. 

Three-body decays are further suppressed with respect to previous formula by another 
power in coupling constant and phase-space reduction, resorting again in the keV range 
of energies. 

From [21] one sees that in most cases single quark decays (SQD) are by far the most 
important decay channels of the bound state, to the order of hundreds of MeV, while as 
discussed above bound state decays are essentially negligible. Moreover a comparison of 
those SQD widths with the results of Table I through eq. (1) shows that for the considered 
mass range of KK there is formation of the bound state. 
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4. Detection 

As wc have previously seen, for large values {R^^ > 300 GeV) SQD is the 
dominant decay channel for a KK bound state, thus leading a to a dominant signature 
consisting of two monochromatic quarks plus missing energy. Following [28] we limit our 
analysis to the region above threshold, i.e. E > 0. The region below threshold, E <(], 
is characterized by peaks in the cross section for values of E equal to binding energies of 
the bound states. The width of those peaks are given by the decay width of the bound 
state, which are at most of the order of the MeV for the SQD and much less, of the order 
of the keV, for other annihilation decay modes, as discussed in sect. 3. 

Because of ISR and beam energy spread, of the order of the GeV for a future linear 
collider, it is unclear whether it could be possible to resolve those peaks of keV magnitude 
with this machine. The only potentially detectable peaks should be the ones belonging 
to a SQD, provided one has a scenario with widths of the order of the MeV. 

The situation above threshold changes drastically with respect to the "naive" Breit- 
Wigner estimate, as is clearly shown in Fig. 1. A few GeV above threshold make for a 
factor of 3 of increase compared to the Born cross- section, allowing a clear distinction 
between the two cases. Assuming an annual integrated luminosity of Lq = 100 fb~^ and 
a scale of the extra-dimension = 400 GeV one finds around 1.2 x 10* events per year 
of two quark decay for a center of mass energy of 10 GeV above threshold (we adopt 
here the scenario for which the branching ratio of SQD is essentially 1). The number of 
events per year loses an order of magnitude at R~^ = 600 GeV, that is about 4 x 10^, as 
could be inferred from Fig. 1. 

For our uiui bound state there are two possible scenarios of decay pattern [24]. The 
first one concerns the iso-singlet Ui^ for which the decay channel into Wi is forbidden 
while that into Zi is heavily suppressed ZiUq^^) ^ sin^ 0i sa 10~^ 10~^ and 

the dominant channel is given by — >■ uo^ji, with — >■ mob7i) ~ 0.98 whose 
signature is a monochromatic quark and missing energy of the KK photon, the latter 
being the LKP [24]. 

For the iso-doublet uij^ the situation is more interesting, as more channels are avail- 
able [24], notably Ui^ — >■ doj^Wi with B{ui^ — )• d^j^Wi) w 0.65 and Ui^ — )■ uo^Zi, with 
B{ui^ — > uoj^Zi) « 0.33 while the branching ratio into 71 is negligible B{uij^ — > 77,0^^71) ~ 
0.02. The decay chain into Wi can follow the scheme: ui^ — )• do^Wi — )• do^£o>^i^ — 
doL^ot^oli with branching ratio given by: 

^ dojoi^oji) « B{ui^ ^ do^Wi)B{Wi ^ loi^i)B{ui ^ i/o7i) « 0.65 i 1 « lO'^ 

and alternatively, the same final state could be reached by the scheme: ui^ do^ Wi — )• 
doj^iii^OL ^ rfot^ot'Ot7i- As compared to the iso-singlet case, the result is a monochro- 
matic quark, a lepton and missing energy in both cases. 

The decay into the Zi channel is ui^ — )• uo^Zi uol^o^i ^01,^0^071, resulting 
in a monochromatic quark, two leptons and missing energy. The branching ratio of the 
above chain is: 

B{ui^ ^ U0M71) « ^ uo^Zi) B{Zi ^ LoLi) B{Li ^ ^71) ~ ^ ^ 1 ~ 5x10"^ 

These leptonic decays of ui have much cleaner signatures than the hadronic ones allowing, 
in principle, for a better detection of the signal. 
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In all cases we emphasize that the observable signal of the bound state production at 
the linear collider would be similar to that of the Born pair production except for the 
absolute value of the cross-section. 

The case of an iso-singlet bound state (or Born pair production of wii?) would give 
rise to the signal e+e~ — )• 2 jets + ^ with cross section: 

(7) cT(e+e" -J> 2jets + aj3,cK x [^{uir wo7i)]^ • 

We note that the ctbick: the iso-singlct ui has to be computed ex-novo and cannot be 
read from the values of Fig. 1 since it refers to the iso-doublet Ui. The singlet and doublet 
have, when including radiative corrections, different masses and the corresponding pair 
production threshold is therefore different. See ref. [34] for details. 

At an e+e~ collider this signal has a standard model background from ZZ production 
with one Z decaying to neutrinos and the other decaying hadronically. Thus for the SM 
background we have: 

(8) £7sM(2jets + 15 =£7zzfbx 0.7x0.2 

In the case of an iso-doublet bound state (or Born pair production oi Ui) the Wi decay 
chain gives the signal: e+e~ 2 jets + 21 + with cross section: 

(9) a{e+e- ^ 2^ + 2£ + = ae^^ [B{uil ^ dojoi^oli)? 

while the Zi decay chain gives rise to the signature e+e~ 2 jets + Ai+ ^ with cross 

sections: 

(10) a{e+e- ^2]+U + ^ = gb^^ [B{uxl ^ wo^447i)]' 

Triple gauge boson production, WWZ, ZZZ at a high energy linear collider has been 
studied in refs. [36, 37]. It has been found that these processes receive a substantial 
enhancement in the higgs mass range 200 GeV < ruH < 600 GeV particularly the ZZZ 
channel. As these processes provide a source of standard model background for our signal 
we estimate them both at a value of nih = 120 GeV and at a value of m/j = 200 GeV for 
which the cross sections are enhanced. Production of WWZ can for instance give rise to 
the signature of 2jcts + 2£ + f via leptonic decay of the W gauge bosons and hadronic 
decay of the Z boson, while the ZZZ production can produce 2jets-|-4£-|-.^via hadronic 
decay of one Z while the others decay leptonically with one of them to a pair of r which 
subsequently decay to ivv {£ ~ e,fi). Estimates of the resulting cross sections are found 
using the CalcHEP [32] and CompHEP [38] software. We have verified agreement with 
previous results given in ref. [37] . We thus estimate within the standard model: 

c7sM(2j + 2^ + 1) « awwz x (0.1)^ x 0.7 

(11) asM(2j + 4£ + 15 « azzz x (0.3)^ x 0.7 x (0.17)^ 

The 2jets + 2£ + ^ channel could be potentially contaminated also from tt pair pro- 
duction cross section which at such high energies is 0(300) fb [39]. Assuming the top 
quarks to decay with probability one to Wb and then the W gauge boson decay via the 
leptonic mode (with B{W — >■ ivi) ~ 0.1) would mimic the signal with a cross section 
o'SM (2jets + 21 + fJ) w 3/6. However in this case we expect 6-tagging of the hadronic 
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R-^ (GeV) 


niui (GeV) 
(iso-singlet) 


2jets + |; 


2jets + 2^ + |; 
ruh = 120 (200) GeV 


2jets + 4£ + |; 
rrih = 120 (200) GeV 


400 


469.0 


81.1 


8.2 (8.1) 


5.6 (5.6) 


600 


703.5 


44.4 


4.1 (4.1) 


3.3 (3.3) 


800 


938.0 


29.0 


2.4 (2.4) 


2.3 (2.3) 


1000 


1172.5 


20.8 


1.6 (1.6) 


1.7 (1.7) 



Table II. - The statistical significance SS as defined in the text corresponding to the annual 

integrated luminosity Lo — 100 fb~^ for the three channels discussed in the text as a function 
of R^^ and ^/s — 2rnui + E, assuming an energy offset of E = 10 GeV from the threshold. In 
the second column we give the values of the ui iso-singlet level-1 KK quarks whose m,asses are 
different form those of the corresponding U\ state from Fig. 1. For the two multilepton channels 
SS has been computed for two values of the Higgs mass ruh = 120 (200) GeV (AR = 20). 



jets. Assuming an efficiency in 6-tagging of 60% wc would get a contribution of 1.2 fb to 
the 2jcts + 2^ + cross-section which has to be added to that in Eq. 11. This has been 
done in the calculation of the statistical significance of tabic II. 

We conclude providing an estimate of the statistical significance 5* 5* = Ns/\/Ns + A^b, 
of the three signals discussed above as related to an integrated luminosity of Lq = 100 
fb~^ {Ng is the number of signal events and A^b is the number of background events). 
These estimates are given in table II. Albeit quite encouraging (especially so the SS of 
the 2jets + ^ wc should bear in mind that the actual observation of these signals might 
be not be so easy from the experimental point of view. Indeed it is quite likely that in a 
framework of a quasi degenerate KK mass spectrum the jets will be typically quite soft 
and therefore difficult to detect. It is therefore customary to concentrate on the much 
cleaner multilepton signatures [24, 40]. An analysis similar to the one given here, but 
with a perspective on signals arising at the Compact Linear Collider (CLIC), regarding 
the (Born) pair-production of level-1 KK- leptons and level-1 KK- quarks is given in [40]. 



5. — Conclusions 

We have considered the formation and decay of a bound state of level-1 quark Kaluza- 
Klein excitation in UED and its consequent detection at a linear e+e" collider. Since 
niKK should be larger than at least 300 GeV we have used a model with a Coulombic 
Being a bound state we have used the Green function technique for the evaluation of 
its formation cross-section in the threshold region, which is more appropriate than the 
standard Breit-Wigner picture as it takes into account the binding energy and the peaks 
of the higher level excitations that coalesce towards the threshold point. The net effect 
is a dramatic increase of the cross-section in the continuum region right of the threshold. 
This multiplicative factor is roug hly 2.6 for R-^ = 400 GeV and drops down to 2.2 at 
= 1000 GeV. The Green function cross-section would allow more than w 10* events 
per year even at = 400 GeV {mu-^ « 478 GeV) for a suitable integrated luminosity 
of the e+e" linear collider {Lq = 100 fb~^). The number of events at R~^ = 1000 GeV 
{mui ~ 1200 GeV) would still be w 10^ at the same integrated luminosity. 

The large difference among the two descriptions of the cross-section should also pos- 
sibly help in the determination of the correct model for such a heavy bound state outside 
the SM. 
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Our analysis of the backgrounds to the final states signals, though very simplified, 
indicates that the multi-lepton channels have a good statistical significance {SS > 2) 
at least up to = 600 ^ 700 GeV, which certainly warrants further detailed and 
dedicated studies of these channels and their backgrounds. The potentially larger (by one 
order of magnitude) statistical significance of the 2j + f channel must be taken however 
with great caution because this signal may be difficult to observe as it is characterized by 
soft jets within the relatively degenerate mass spectrum of the extra-dimensional model. 
Further detailed studies are also needed for this channel. 

* * * 

N. F. was supported by the Fondazione Cassa di Risparmio di Spoleto. 
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